-->[f,xopt,gradopt,work]=optim(list(purequad,Qmat),[1 1],imp=9)
 norm of projected gradient lower than   0.0000000D+00
 
 work  =
 
 
         column 1 to 6
 
    6.0000008    2.950D-09    78.    1.755-163    3.043-166    1.053-162  
 
         column  7 to 12
 
    2.373-164  - 1.755-163  - 3.043-166  - 1.053-162  - 2.373-164  - 2.297-170  
 
         column 13 to 15
 
  - 3.982-173  - 1.378-169  - 3.106-171  
 gradopt  =
 
   1.0-169 *
 
  - 1.3781362  - 0.0310601  
 xopt  =
 
   1.0-170 *
 
  - 2.2968937  - 0.0039821  
 f  =
 
    10.  
 
-->diary(0)